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Abstract 

^' 

The irreversibility of the dynamics of the conservative systems on example of hard disks 

and potentially of interacting elements is investigated in terms of laws of classical mechanics. 

The equation of the motion of interacting systems and the formula, which expresses the 

entropy through the generalized forces, are obtained. The explanation of irreversibility 

mechanism is submitted. The intrinsic link between thermodynamics and classical mechanics 

£3 was analyzed. 

a 

O ' Introduction. 

Irreversibility determines the contents of the second law of thermodynamics in fundamental physics. 

According to this law there is function S named entropy, which only grows for isolated systems, achieving 

> ■ 
^\ • a maximum in an equilibrium state. But the potentiality of forces for all four known fundamental 

—*. ' interactions of elementary particles is a cause of reversibility of the Newton equation [1] and therefore 

O, 

parts contradicting each other: there are reversible classical mechanics and irreversible thermodynamics 



o 
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The first attempt to resolve these contradictions belongs to Boltzmann. From the 7J-theorem it 
follows that the many-body systems should equilibrates. But at reception of these results Boltzmann used 
probabilistic principles. Therefore the contradiction between classical mechanics and statistical physics 
g i has not been overcome. Great importance of this problem for all physics explains its big popularity 
O ■ among the famous physicists of the world. The history of its solution is very extensive and at times 



dramatic. Therefore let us reference only several works, which precisely enough describe a modern 
condition of a problem of irreversibility [2-4] . 

For overcoming of the irreversibility problem I. Prigogine has suggested to try to create the expanded 
formalism of open systems within the framework of classical mechanics laws [4]. Indeed, by analyzing 
dynamics properties of a hard disks system it turned out that such formalism appears in a process of the 
solution of a problem of irreversibility [5-9]. For example, as a result of this analyzing the generalized 



Liouville equation and equation of the motion for the systems was obtained. Here on the background 
of results [5-9] the explanations of the mechanism of irreversibility, interrelation of classical mechanics 
and thermodynamics are offered. All this determines a way of construction of the expanded formalism 
of open systems within the framework of classical mechanics laws. 

Our investigation is based on the next method. The conservative system of interacting elements is 
prepare in the noncquilibrium way. This system is splitting on subsystems so that they could be accepted 
as in equilibrium. The subsystems dynamics under condition of their interactions is analyzed. 

First of all we investigate the dynamics of the simplest hard-disks systems. For those we find that 
so-called generalized forces between selected subsystems should depend from velocities of those. This 
leads to the key idea about impossibility of solution of the problem of irreversibility in terms of canonical 
Hamilton or Liouville equations. So, we have built the generalized Hamilton and Liouville equations. 
Then we deduce how the irreversible dynamics follows from the generalized Liouville equation. 

After that we analyze the dynamics of Newtonian systems and deduce the equations of the motion 
for the underlying subsystems. On the base of these equations we answer the question: why and how 
the velocity dependence for the generalized forces between the subsystems appears for the case of the 
potential-type forces between the elements. Next point: the intrinsic link between thermodynamics 
and classical mechanics was analyzed. The general formula for the entropy of the system in terms of 
generalized forces is obtained which corresponds to irreversible transformation of the interaction energy 
of subsystems into the internal energy by means of the generalized forces work. 

Irreversibility in a hard- disk system. 

We used an approach of pair interaction of disks. Their motion equations are deduced with help of 
the matrix of pair collisions. As complex plane this matrix is given [5]: 

(a —ib \ 
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where a = dkj exp(ii?/jj); b = (3 exp(i$kj)', dkj = cos'dkj', (3 = sinfikj] i- an imaginary unit; k- 

and j- numbers of colliding disks; dkj- the impact parameter (IP), determined by distance between 
centers of colliding disks in the Cartesian plane system of coordinates with axes of x and y, in which 
the /c-disk swoops on the lying the j- disk along the x - axis. The scattering angle "&kj varies from 
to 7r. In consequence of collision the transformation of disks velocities can be presented in such form: 
Vki = SkjVkj ( a )j where Vj~ and V k + - - are bivectors of velocities of k and j - disks before (— ), and after 
(+) collisions, correspondingly; Vkj={Vk,Vj}, Vj = Vj x + iVj y - are complex velocities of the incident 
disk and the disk - target with corresponding components to the x- and y- axes. The collisions are 
considered to be central, and friction is neglected. Masses and diameters of disks-" d" are accepted to 
be equal to 1. Boundary conditions are given as either periodical or in form of hard walls. From (a) we 



can obtain equations for the change of velocities of colliding disks [5-7]: 

Here, A k j = Vk — Vj - are relative velocities, 5V k = V k + — V k ~ , and 5Vj = V + — V~ - are changes of disks 

velocities in consequence of collisions, ip k j = if3e ki . 

That is, eq. (1) can be presented in the differential form [5]: 

V k = $ kj 5(ip kj (t))A kj (2) 

t 
where i/'fej = [1 — Kfej|]/I^fej|i ^(V'fcj)-delta function; l k j(t) — zjL + J A k jdt - are distances between 

o 
centers of colliding disks; z° = z k — Zj, z® and Zj - are initial values of disks coordinates; Q k j = 

i(l kj A kj )/(\l kj \\A kj \). 

The equation (2) determines a redistribution of kinetic energy between the colliding disks. It is not a 

Newtonian equation because the forces between the colliding disks depend from their relative velocities. 

So, we get the generalized Hamilton equation to be applied for studying the subsystems dynamics [6] : 

^--A + iT 0) 

^ - 14 (4) 
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These are the general Hamilton equations for the selected p-subsystem. The external forces, which acted 
on p-subsystem, presented in a right-hand side an eq. (3) . 

Using eq.(3,4), we can find the Liouville equation for p-subsystem. For this purpose, let's to take a 
generalized current vector - J p = {r'k,pk) of the p-subsystem in a phase space. From equations (3,4), we 
find [6,7] : 

d ^=->>iM <5> 

Eq.(5) is a Liouville equation for p-subsystem. It has a formal solution: 
/„ = cons* -exph/d: ^F^dt}. 

fe=l 

The equation (5) is obtained from the common reasons. It is suitable for any interaction forces of 
subsystems. Thus, the eq.(5) is applicable to analyze any open nonequilibrium systems. In particularly, 
it can be used for explanation of irreversibility. 

T 

The right side of eq.(5), f p Y^ w~^k^ IS * nc integral of collisions. This integral can be obtained from 

fe=i Pk 
the motion equations of the systems element. For example, for a hard disks system it can be found with 

the help of the eq. (2). 

Let's consider the important interrelation between descriptions of dynamics of separate subsystems 

R T 

and dynamics of system as a whole. As the expression, ^ ^ F£ = 0, is carried out, the next equation 

p=i fc=i 



for the full system Lagrangian, L,r, will have a place: 

d 8Lr OLr 



(6) 



dt dV k dr k 

and the appropriate Liouville equation: -4jf- + Vk-^ 2 - + Pk~Q = 0. The function, /#, corresponds to 

R 
the full system. The full system is conservative. Therefore, we have: E divJ p = 0. This expression is 

P =i 

R R R ~ 1 R R 

equivalent to the next equality: ^(E ln /y) = ft( ln U. f P ) = (U. f P ) ft(U f P ) = °- So > U fp = 

p— l p— l p— l p— l p— l 

R R 

const. In an equilibrium state we have J\ fp — Ir- Because the equality E Fp — i s fulfilled during all 

p— l p— l 

time, we have that equality, J~J f p = / fl , is a motion integral. It is in agreement with Liouville theorem 

p=i 
about conservation of phase space [10]. So, only in two cases the Liouville equation for the whole 

system is in agreement with the general Liouville equation for selected subsystems: if the condition 

t T " ' T 

J (E w~Fk)dt — ► const (c) is satisfied when t — > oo, or when, (E w~^k)^ IS a periodic function 

o fc=i Pfc fe=i Pk 

of time. The first case corresponds to the irreversible dynamics, and the second case corresponds to 

reversible dynamics. Because the generalized forces for a hard-disk system depended from velocities, the 

irreversible dynamics is possible. 

Dynamics of strongly rarefied systems of potentially interacting elements is also described by the eq. 
(2). Therefore for those systems the irreversibility is possible as well [7,8]. 

The dependence of generalized force from velocities - it is a necessary condition for the irreversibility 
was really to be occurred. Therefore the question about irreversibility for Newtonian systems is reduced 
to that about the velocity dependence of forces between subsystems. 

For hard disks and strongly rarefied systems of potentially interacting elements the presence of 
irreversibility is predetermined by equation of motion (2) . In these systems the interaction forces between 
the elements is depending from velocities. Therefore it is clear that the generalized forces will depend 
from velocities also. But forces between elements for Newtonian systems arc potcntional. Therefore it is 
necessary to answer the question: how velocity dependence of generalized force is appearing when forces 
between the elements are independent from velocities. For this purpose in the next part we will obtain 
of the motion equation for subsystems. 

The evolutionary equation of the motion for subsystem. 

Let us to analyze Newtonian systems. We take a system with energy: En = TV + Un = const, 

N 

where TV = \ E v i 2 ~' IS a kinetic energy; f7jv(^ij) - is a potential energy; r^ = r« — rj - is a distance 
between i and j elements; N - is a number of elements. Masses of elements are accepted to 1. 
The equation of the motion for elements is: 

N „ 

* = - E ^ (7) 

It is reversible Newton equation. The irreversibility in Newtonian systems can be compatible with 
reversibility of the Newton equation in the case if exists the parameter determining irreversibility, con- 



corning which Newton equation is invariant. Such parameter is a velocity of the center of mass of 
subsystem. The Newton equation for subsystem elements is invariant in respect to this velocity. But 
an exchange of energy between subsystems is not invariant to this velocity. These energy exchanges are 
determined by the work of the generalized forces, which are proportional to the speed of change of the 
velocity of the center of mass of subsystem. 

Easiest explanation of the fact, why the Newton equation does not determine the generalized forces, 
can be made on the example of system consisting of two balls connected by a spring. Shortly, in the 
laboratory system of coordinate the energy of each ball simultaneously transforms in kinetic energy of 
other ball and in the potential spring energy. But the Newton equation describes only transformation 
of kinetic energy into the potential energy. Therefore to solve this task with the help of the Newton 
equation, it is necessary to exclude energy of the system motion as the whole. For two bodies it is 
possible to do by transition to the coordinates system of the center of mass. In this coordinate system 
the change of kinetic energy of a ball is equally to the change (but with opposite value) of the potential 
energy. But for three and more bodies it is impossible to do. Therefore N > 3 bodies system in general 
case cannot be reduced to the independent equations of Newton [4,8-10]. 

Therefore, for the description of evolution of nonequilibrium system it is necessary to obtain the 
equation determining generalized forces or the speed of change of the velocity of the center of mass 
of system, in results of its interaction with other subsystems. For this purpose we'll do the following. 
In the some system of coordinates we shall present kinetic energy of system as energy of subsystem 
motion as the whole - T^ r , kinetic energy of motion of its elements concerning the center of mass -T^ s , 
a potential energy of subsystems elements- U l ^ s and interaction energy with other subsystems. The 
energy E l £ s = T"^ s + U]^ s we'll call a binding energy. At absence of external forces the energy Tjy 
and E™ 8 are the motion integrals. Their sum is a full energy of the system. We shall assume, that the 
system is prepared by nonequilibrium way so, that it can be divided into two equilibrium subsystems. 
Then it is possible to generalize results of their analysis on any nonequilibrium systems. 

The following equations for exchange of energy between subsystems have a place [9] : 

«ife + t EKll + ££,}-£ XXg£L (8) 

j=i+l i=l l ° j K =li L =l % ^3K 

K K ~ x arr k l „ 

KVkVk+Y. EW? + & B -SE^ W 

j=i+li=L+l %3 j K = li L = l lL3K 

Here we take LVl + KVk = 0, Vl and Vk is a velocities of the center of mass of L and K subsystems; L 
and K are number of elements in subsystems; i>ij -is a relative velocities % and j elements; L + K = N . 
Masses of elements are accepted to 1. Double indexes are accepted for a designation, to what subsystems 
these elements belong. 

The left side in eq.(8,9) determines exchanges of the energy between subsystems. The first members 
set change kinetic energy of motion of subsystems, T tr as the whole. The second members describe 



transformation of binding energy of the subsystems, E ma . 

The right side in the equations eq. (8,9) determines potential energy of interaction of subsystems. 
They cause transformation of kinetic energy of the motion of subsystems, T tr , to their binding energy. 

In case L = K the equation of the motion for one subsystem can be obtained from eq.(8,9) [9]: 



l l-\ . „ rr K L 

Vl = -v7 L E E^if + d>-EEs- do) 






The equation (10) is determining the generalized force. This force is depending not only from 
coordinates but from velocities also. Presence of the dependence from the velocities is a necessary 
condition of irreversibility. Therefore we will call this equation as evolutionary equation. With the help 
of eq.(2) it is not difficult to show that equation (10) is correct for hard-disks system. 

In a limit N 3> 1 for enough fine-grained splitting system, the equation (10) determines a field of 
generalized forces between the subsystems, caused by their relative motion. Irreversibility is connected 
with decrease of the generalized forces in a result of transition of kinetic energy of motion of a subsystem 
in the binding energy of subsystems. For any splitting of the equilibrium system the relative velocities 
of subsystems and the energy flow between them are equal to zero [9, 11]. 

Let us explain, why the generalized forces depend from velocities of the elements. Inherently, this 
force determines the change of velocity of the center of mass. But this velocity is determined by distri- 
bution of elements velocities. Therefore the interaction force of subsystems should depend not only from 
the forces of potential interaction of elements, but also from distribution of their velocity. 

Laws of conservation of energy and momentum do not forbid for equilibrium system of spontaneous 
occurrence of an inequality T tr > for some of subsystems. This inequality means infringement of 
equilibrium [10, 11] and transition from canonical Hamilton formalism systems description to the de- 
scription of system inside the framework of the generalized equations. It is connected with appearance 
of the generalized forces between subsystems. But canonical Liouville equation cannot be transformed 
by spontaneous way to the generalized type. It is impossible because of conservation of the canonical 
Hamilton's formalism. From the physical point of view it is impossible because of stability of an equilib- 
rium state for mixed systems [8] . Such stability is caused by aspiration to zero of the generalized forces 
arising at a deviation of system from equilibrium [8, 11, and 14]. Hence, the system, having come to 
equilibrium, never leaves this condition. 

Thermodynamics and classical mechanics. 

Let's consider interrelation of classical mechanics with thermodynamics. Existence two invariant of 
the motion- E lns and T tr , and also character of their transformation at interaction of the subsystems, 
determined by the equation (10), allow to catch deep analogy between the equations (8-10) and the basic 
equation of thermodynamics [14]: 

dE = dQ-PdY (11) 



Here, according to the common terminology [14], E is internal energy of a subsystem; Q is a thermal 
energy; P is a pressure; Y is a volume. 

The change of energy of the selected subsystem is caused by work of external forces. Therefore to 
change of full energy of a subsystem corresponds to dE. 

The change of kinetic energy of motion of a subsystem as the whole, dT tr , corresponds to the member 
PdY. Really, dT tr = VdV = VVdt = Vdr = PdY 

Let's us determine, what the member in the eq.(ll) corresponds to the change of the binding energy 
of a subsystem. From the virial theorem [8] it follows that if potential energy is homogeneous function 

_ . — ins - ins 

of the second degree from all radiuses-vectors, then E ms = 2T = 2U . Feature means averaging on 
time. Earlier we have obtained, that the binding energy, E ms increases due to the energy, T tr . But 
opposite process is impossible. Therefore the member Q from the eq.(ll) is conforming to the energy 

rpins 

Let the subsystem will consist of N m elements. Then average energy of each element, Tq™ s = E ms /N m . 
Let the binding energy increases on the value dQ. In connection with virial theorem, keeping members 
of the first order, we shall have: dQ ~TQ U8 [dE ms /Tl ns ] — T™ s [dv/vQ\ 1 where «o is average velocity 
of an element, and dv its change. As the condition the subsystems equilibrium is accepted, we have 
dv/vo ~ dT m /T m . Here T m is a phase volume of a subsystem, and dT m is increasing T m due to in- 
creasing of the subsystem energy on the value, dQ. By keeping members of the first order we'll have: 
dQ w KE ins dT m /T m = nE ms d\nT m , where T ms = nE ms . By definition d\nT m = dS lns , where S ms is 
entropy [10,14]. So, we have near equilibrium: dQ « T^ ns dS ms = nE ms dS ms . 



Let's consider the connection of the field of forces of interaction of subsystems with entropy. As the 
sum of all forces of interaction of subsystems at any moment is equal to zero, the condition, ^ T^ = 



m— 1 



have a place. I.e. energy T is redistributed between subsystems and goes on increasing of the binding 
energy. It is equivalent to entropy increasing. Relative velocities of subsystems go to zero, and the system 
equilibrates. Near to equilibrium the right side of the eq.(5) is equal to zero. So, entropy deviation from 
equilibrium state can be determined with the help of following formula: 



A5 = E{JS7E/E^ m ^} ( 12 ) 

1=1 k=l J s 



Here E m ' is a kinetic energy of subsystem; mi is a number elements in subsystem -" V ; R- is a number 
of subsystems; s- is external disks which collided with internal k disk along trajectory. The integral is 
determining the work of the force F™ s l during equilibrating. 

The eq.(12) connects the generalized force with entropy. Similar connection is established by the 
formula for KS - entropy in which entropy is expressed through Lyapunov's exponents [3]. Inherently the 

R 

eq.(12) corresponds to the formula for entropy: S = J2 Si{El ns + T/ r ), (see, [11]). Really, if E\ ns » T/ r , 

i=i 

R 

then dS =J2 -§r& dT i r that corresponds to eq.(12). 
l=i l 
Summarizing. 



Thus, the dynamics of the open and noncquilibrium systems is determined by the generalized Li- 
ouville equation and the evolutionary equation of the subsystems motion. The necessary condition of 
irreversibility follows from generalized Liouville equation. This condition is dependence of the general- 
ized forces from relative velocities of subsystems of any non-equilibrium system. So, due to the work of 
the generalized forces, the energy of relative motion of subsystems is transformed to their binding energy 
by irreversibly way and goes on entropy increasing. But the binding energy cannot be transformed to 
the energy of relative motion of subsystems. This is essence of the mechanism of equilibration [5-9, 11]. 

Though dynamics of hard disks can be described without use of potential energy, the essence of the 
mechanism of equilibration for them is identical to the mechanism for Newtonian systems. Difference 
only in that that in Newtonian systems the energy of relative motion of subsystems can be transformed 
to their binding energy due to the potential interaction of subsystems. 

In equilibrium we have TjJ'=Tj[=0 and Vl = Vk = 0, therefore the work of interaction force of each 
subsystem, and also energy T tr are equal to zero. It is equivalent to conditions: E™ 8 = E l g s = 0. Hence, 
the closeness of system to equilibrium means the closeness to zero of the generalized forces between sub- 
systems. The criterion of equilibrium is applicability of the canonical Hamilton and Liouville equations 
for the description of dynamics of subsystems for any splitting of full system. Only in equilibrium the 
dynamics is reversible for both the subsystems and for the system as a whole. The dynamics is reversible 
also for the weak fluctuations because the entropy has a second infinitesimal order in respect to deviations 
of a system from equilibrium [9, 11, 13]. 

Criterion for the sufficiently fine-grained splitting of system is the closeness of relation T tr /E ms to 
zero. It is obvious that the number of elements in a subsystem should be large enough. Otherwise 
fluctuations become commensurable with errors of a field of generalized forces [11, 13]. 

Splitting the system onto the equilibrium subsystems is necessary reception for studying the dynamics 
laws of the many-body systems. It allows setting a field of the generalized forces and to describe it's 
with the help of the evolutionary equation (7). In particular, only due to this reception it is possible to 
find out, that their work is goes by irreversible way on increase of the binding energy of a subsystem. 

The basic results received here do not contradict results of the analysis of systems, which was 
obtained in earlier works by similar method: splitting into the subsystems and studying its dynamics 
near equilibrium. So, in [11] the macroscopical motion (it is similar to the subsystems motion in our 
case) was considered. There by using of the method of multipliers Lagrange it has been strictly proved 
that equilibrium is possible only when T tr — 0. Also the validity of the equation for entropy (12) is 
follows from it. 

Why till now it was impossible to solve a problem of irreversibility? From our point of view the 
explanation is consists in the following. As a rule it was accepted without doubt that the full description 
of evolution of noncquilibrium systems could be realized on the basis of the Newton equation and 
canonical Liouville equation. But we find, that the canonical Liouville and Newton equations are not 



applicable for describing work of the generalized forces because in nonequilibrium systems they are 
dependent from velocities. 

Indeed, let us try to determine the generalized force with the help of Newton eq.(7) by summation 
of all forces acting on elements of a subsystem. Then we obtain the forces independent from velocities 
due to invariancy of the Newton equation concerning Galileo transformation . So, the Newton equation, 
does not describe a field of the generalized forces. 

Let's to emphasize, that the substantiation of irreversibility mechanism and connection classical 
mechanics and thermodynamics offered here, follows only from laws of conservation of energy and a 
momentum. 

The evolutionary equation of the motion subsystem along with the generalized Hamilton and Liouvillc 
equations given above, and also the established relations between classical mechanics and thermodynam- 
ics, can be put into the base of theoretical description of evolution of the open systems. 
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